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We consider a simplified constituent model for relativistic strong-interaction decays of hybrid 
mesons. The model is constructed using rules of renormalization group procedure for effective parti- 
cles in light-front quantum field theory, which enables us to introduce low-energy phenomenological 
parameters. Boost covariance is kinematical and special relativity constraints are reduced to the 
requirements of rotational symmetry. For a hybrid meson decaying into two mesons through disso- 
ciation of a constituent gluon into a quark-anti-quark pair, the simplified constituent model leads to 
a rotationally symmetric decay amplitude if the hybrid meson state is made of a constituent gluon 
and a quark-anti-quark pair of size several times smaller than the distance between the gluon and 
the pair, as if the pair originated from one gluon in a gluonium state in the same effective theory. 

PACS numbers: 13.25.Jx, 11.80.-m, 13.90.-H 



I. INTRODUCTION 

Hybrid mesons, in particular those with exotic quan- 
tum numbers, will play important role in the phe- 
nomenology of confinement and studies of low energy 
gluons. For example, decay patterns of gluonic excita- 
tions in hybrid mesons are expected to be quite different 
from those that characterize ordinary meson resonances. 
It is predicted that exotic meson decays are dominated 
by the so-called "S+P" channels that contain pairs of 
mesons, one with the valence quark-anti-quark pair in 
the S-wave relative motion, and the other with a pair in 
a P-wave, e.g., nbi, nfi, or ixa\ Q. In contrast, non- 
exotic meson decays are dominated by "S+S" channels, 
such as 7T7T, pp, KK, etc. These predictions are based 
on a simple, non-relativistic constituent-quark picture. 
However, the lightest exotic mesons are expected to be 
as heavy as 1.6 — 2 GeV and their decays will produce 
mesons that move relativistically. Therefore, this article 
addresses the issue of constraints that special relativity 
imposes on the constituent picture of such decays. The 
other basic requirement satisfied by the relativistic ap- 
proach developed here is that in principle it is derivable 
from quantum field theory in agreement with rules of 
renormalization group for low energy parameters in the 
constituent Hamiltonians. 

In order to analyze data from experiments where pho- 
tons are absorbed on a hadronic target and produce in- 
termediate hadronic states with masses on the order of 
1.5-2.5 GeV, which subsequently decay into mesons 7r, p, 
or b, one must construct a viable picture of strongly inter- 
acting bound states. These states, especially 7r-mesons, 
move with velocities comparable to the speed of light, 
and the whole photo-production and decay process can- 
not be described using wave functions for hadrons at rest, 
or in slow motion. At the same time, the principles of 
quantum field theory have to be obeyed precisely if one 
demands a clear connection between a phenomenologi- 



cal constituent picture of the process and QCD. Lessons 
from QED are of great value but QCD interactions are 
much stronger and pose additional problems with non- 
perturbative and relativistic effects, which QED is not 
telling us how to solve beyond Feynman's perturbation 
theory. 

So far, theoretical insights into the exotic mesons struc- 
ture are provided by lattice simulations , and Coulomb- 
gauge QCD ||, in addition to model studies || |j. It 
follows from these studies that the lowest mass and spin 
exotic mesons are most probably dominated by a con- 
stituent gluon coupled to a constituent quark-anti-quark 
pair. In a more complete picture, one would like the 
exotic candidate states to be obtainable, at least in prin- 
ciple, from a precisely defined relativistic Hamiltonian 
for QCD. The decay process should be described by the 
same Hamiltonian. To satisfy these requirements, one 
would have to derive the effective constituent particles 
from quantum field theory. The derivation would have to 
include description of the constituents' spin and angular 
momenta, and how rotating and boosting of bound states 
to high velocities is accomplished. Such comprehensive 
relativistic constituent picture is not yet established in 
QCD and hadronic constituents are far from being un- 
derstood as well as the non-relativistic constituents of 
atoms are in QED. However, there exist reasonable the- 
oretical grounds for hope that the required relativistic 
description of effective quarks and gluons can be found 
using renormalized light-front (LF) QCD Hamiltonians 
in A + = gauge. The reasons for the hope go be- 
yond the pioneering works of Lepage and Brodsky Q and 
Wilson et al. [||, and include results that suggest that 
the effective constituents are derivable in a perturbative 
renormalization group calculus for the corresponding cre- 
ation and annihilation operators |8j. The same calculus 
is also suggested to be able to produce all generators 
of the Poincare algebra in the form ready to use in the 
effective-particle Fock-space basis H. If these genera- 
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tors were available, the desired construction of hadrons, 
including exotic hybrid states and fast moving mesons, 
could be carried out precisely enough to make specific 
predictions for partial waves branching ratios in decay 
amplitudes, and the QCD constituent theory of hybrids 
could be put to definite tests. 

The LF approach is a potential alternative to the stan- 
dard form of Hamiltonian dynamics, where states evolve 
in usual time. There, rotational symmetry is kinemat- 
ical, but one needs to worry about boosting the con- 
stituent wave functions for the fast moving bound states 
and their constituents fill . Nevertheless, much progress 
has been achieved in building the constituent representa- 
tion in the usual Coulomb-gauge quantization in the rest 
frame || [l2j . The connection between bare field quanta 
and quasi-particles, corresponding to a constituent basis, 
was obtained and related to spontaneous chiral symme- 
try breaking. Furthermore, in the gluonic sector, the 
Coulomb-gauge canonical quantization reproduces the 
the lattice gauge theory results for the spectrum of the 
low-lying gluon modes, e.g., glueballs, and the excited 
heavy-quark adiabatic potential Still, the Lorentz 
boost covariance of observables in the instant formula- 
tion in Minkowski space is not yet established, and poses 
problems analogous to rotational symmetry issue in the 
boost invariant LF scheme. Both are major theoretical 
challenges. Here, we will concentrate on the rotational 
symmetry in the LF description of hadrons in terms of 
their constituents. 

Note, however, that the LF QCD path is not suffi- 
ciently developed for immediate identification of all im- 
portant features of the effective constituent dynamics 
in hadronic processes. The theory requires initial guid- 
ance from phenomenological analysis to begin scanning 
through and cataloging of the infinite set of potentially 
relevant operators that are allowed to appear in the ef- 
fective QCD Hamiltonian. In view of the magnitude of 
the task, before one dives into the full complexity of LF 
QCD, one should first find out if rotational symmetry of 
observables has any chance to be obtained in the boost- 
invariant picture of hadronic states made of a small num- 
ber of effective-particle Fock-components. The common 
belief is that a physically adequate relativistic scheme 
cannot be achieved using a small number of degrees of 
freedom. The structural integrity requirement for the 
effective-particle LF-scheme thus implies that the rota- 
tional symmetry condition should be possible to satisfy 
with some natural easiness, visible already in the most el- 
ementary versions of the scheme. Therefore, the question 
of how hard it is to build a rotationally symmetric model 
of a hybrid decay amplitude with a minimal number of 
constituents needs to be answered before a search for ap- 
proximate solutions to quantum field theory in terms of 
effective constituents is pursued on a larger scale. Here, 
we report results of the feasibility test for this scheme, 
which were obtained with all possible simplifications one 
could afford preserving only the most basic and universal 
constraints of special relativity, quantum field theory, and 



renormalization group for effective particles. All details 
beyond the core features are removed for transparency of 
the bottom-line structure. 

In principle, the dynamics of hadronic constituents is 
determined by the effective-particle Hamiltonian. Once 
such a Hamiltonian is constructed, one can attempt to 
solve the corresponding Schrodinger equation and calcu- 
late the hybrid decay in terms of the effective-particle 
interactions and states. The basic study reported here 
is carried out instead by simplifying the initial canon- 
ical Hamiltonian theory to the ultimate minimum, de- 
riving the effective theory in a crudest lowest-order ap- 
proximation, and using model wave functions typical for 
phenomenology instead of self-consistent solutions of the 
currently unknown effective theory. Thus, in practice, 
we calculate a transition amplitude for a hand-made 
one three-effective-particle bound state to decay into two 
hand-made two-effective-particle bound states, and we 
check if the result has any chance to respect rotational 
symmetry. 

This paper is organized as follows. Section II describes 
the key elements of the constituent scheme used here. It 
is assumed that the interaction responsible for the hy- 
brid decay comes from an effective Hamiltonian obtained 
from a local quantum field theory via renormalization 
group procedure for effective particles. Section III de- 



scribes details of the model states of hybrids and mesons, 
and the transition amplitude for the hybrid decay into 
the two mesons. Wave functions of the incoming and 
outgoing states are written down on the basis of well- 
known phenomenology, but without solving the effective 
Schrodinger equation. Nevertheless, we allow a whole set 
of free parameters in them in order to cover a large range 
of possible structures and find out if it is possible (and if 
so, then for what choices of the parameters) to obtain ro- 
tationally symmetric results. Numerical results from our 
review of the angle dependence of the decay amplitude 
are discussed in Section IV. Section |y| draws conclusions 
that concern phenomenological studies of hybrid produc- 
tion and decay processes and can be directly useful in 
comparing data from future experiments with a theory 
of effective particles in QCD. 



II. 



KEY ELEMENTS 



In order to remove the basic issues of spin, gauge in- 
variance, current conservation, and effective mass of glu- 
ons, from the preliminary constituent study, and to focus 
on general quantum mechanical and relativistic features 
of the effective particle dynamics, one can assume that 
gluons are scalars. In order to do better than that at 
the current level of understanding of the effective gluon 
dynamics, one would have to clearly answer questions 
concerning Lorentz transformation properties of massive 
gluons with spin one and only two, instead of three po- 
larization states. Since the answer is lacking, the issue is 
avoided by the arbitrary simplification. 
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With the gluon spin issue removed, the scalar-iso- 
scalar bare gluon can still be coupled to an iso-doublet 
of bare quarks like in a Yukawa theory. Thus, the bare 
"gluon" has only color as in QCD, while the coupling to 
bare quarks is given by the term 

H mt =gM a t a 4> , (1) 

where the gluon field does not have a vector index. In 
this model, the first order solution of renormalization 
group equations for the effective gluon coupling to ef- 
fective quarks produces the term 

gMxA a x t a ^x , (2) 

where A denotes the renormalization group parameter. It 
has dimension of mass. Its scale can be intuitively associ- 
ated with the transverse size of the effective particles and 
fx denotes the vertex form factor. The width of the form 
factor in momentum space is given by A. The effective 
particles that can be considered appropriate for the low- 
est Fock-sector picture of hadrons are expected to emerge 
when A is lowered by solving similarity renormalization 
group equations for Hamiltonians down to the scale of 
hadronic masses, possibly to the values on the order of 1 
to several GeV, when one aims at description of a hybrid 
meson in terms of a quark-anti-quark pair and a gluon. 

The factor g stands for the vertex strength function 
v\ that at some value of the momentum transfer defines 
a running coupling constant in the Hamiltonian de- 
noted by g\. In the first order solution of renormalization 
group equations, v\ equals g. In higher orders, this cou- 
pling would be replaced by g\ and the latter could be 
quite large in QCD when A is of the order of hadronic 
masses. However, the large coupling is not sufficient to 
make the effective interaction strong enough to wash out 
the constituent picture of the individual hadrons, as it 
would happen in a local canonical theory with equally 
large coupling constants. The reason is that the vertex 
form factor cuts off interactions with momentum trans- 
fers larger than A, and they are effectively much weaker, 
like in nuclear physics. Therefore, the bound state pic- 
ture with a small number of massive constituents may 
be still valid when A is small. Eq.(j|) provides also a 
potentially valid structure for at least a part of the inter- 
action responsible for the decay of a model constituent 
"gluon" into the effective quark and anti-quark pair. The 
"gluon" decay factor in the model hybrid decay ampli- 
tude can thus be assumed to be given by the product of 
the form 

fl / A (12)u 1 wa , (3) 

times the corresponding color factor xJ 1 i°Xc2- One may 
simplify the constituent picture even farther for the pur- 
pose of this study, and replace fermionic quarks by 
scalars, too. The required interaction, analogous to (Q), 
would be 

H int =g?A a t a Z , (4) 



where £ is the colored scalar quark field .and the interac- 
tion factor in the decay, analogous to (||), would reduce 
to 

S/a(12) , (5) 

the color factor being the same as for fermions. Both 
cases of spin- 1/2 and spin-0 constituents are described 
in this paper to show that the features they display are 
general and independent of spin. This is one of the rea- 
sons why one can expect them to occur also in the future 
QCD calculations. 

Even a modest attempt at a self-consistent dynamical 
solution of a relativistic field-theoretic model with cou- 
pled three-body and two-body bound states is a large 
project. Our goal here is not to carry out such project 
but to find out if the proposed scheme has a chance to 
produce rotationally symmetric decay amplitude. More- 
over, even if a constituent picture may satisfy the special 
relativity constraints, the correct effective QCD theory 
is not known yet. In order to see how strong the rela- 
tivity constraints are, the simplest thing to do is to as- 
sume that the outgoing mesons and decaying hybrid are 
both scalars and iso-scalars, and to write down candi- 
dates for their wave functions in the effective constituent 
basis with a priori unknown parameters that are con- 
strained only by reasonable guesses about their natural 
sizes. If the rotational symmetry constraints cannot be 
satisfied within this class of models, it is unlikely that a 
simple constituent model can ever be built in more com- 
plex cases with additional dynamical constraints. Thus, 
for example, we assume here only that a meson has a 
Gaussian momentum-space relative-motion wave func- 
tion with width (3. Its spatial radius is of the order of 
1 / (5 (if (3 is not too large in comparison to the constituent 
masses and the meson mass since relativistic effects may 
change connection between (3 and the radius making it 
less constraining |]l3[ ). Thus, in the first approximation, 
a meson of size 0.5 fm is expected to have j3 ~ 0.4 GeV. 
In the theory of scalar effective constituents, no further 
considerations are required for constructing plausible test 
candidates for the relative-motion wave functions, and 
none is incorporated here. 

In the more complex case of fermions, one has to decide 
in addition how to combine spins in the wave functions for 
effective constituents. Solutions of effective fermion the- 
ory are also not known, and one is currently in the similar 
position here as in trying to guess hadronic wave func- 
tions in phenomenological studies in the nonrelativistic 
quark model . The simplest standard phenomenolog- 
ical procedure is to assume that the spin wave functions 
for the lowest effective Fock components are the same as 
for multiplets of Poincare group representations for free 
particles, i.e., products of free Dirac spinors in the case of 
fermions. For scalar (colorless) hadrons made of qq pairs, 
one would use spin factors of the form u\Vi (the color 
factor would be X c iXc2)- The spin amplitudes would be 
multiplied by the relative momentum-dependent factors 
that would be of similar form as in the scalar case. 
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The relative-motion scalar wave-function factors rep- 
resent binding and allow for shifting the total free energy 
of the intermediate constituent particles away from that 
of the energy shell value for the hybrid itself and the 
outgoing meson products. This shift in energy is asso- 
ciated with conservation of the total three-momentum 
of all particles. In the LF Hamiltonian dynamics, how- 
ever, the conserved momenta are not the same as in the 
dynamics developed in Minkowski's time. Instead, the 
conserved momentum components (in conventional no- 
tation) are p + = p° + p 3 and p- 1 , where _L= 1, or 2. The 
role of energy is played by p~ , which for a free particle 
of mass m would be p~ = (p^ 2 + m 2 )/p + . 

The key Lorentz invariance test, associated with the 
complicated dynamics of rotations in the LF coordinates, 
is reduced to the question if the off-shell energy shift in 
the direction of p~ does, or does not introduce significant 
artificial angular asymmetries in the model hybrid decay 
amplitude. For example, can one avoid a sizable depen- 
dence of the amplitude on the angles at which the scalar 
mesons fly away from the scalar hybrid decay event in 
the hybrid center of mass system of reference? If such 
angle dependence would unavoidably develop in the sim- 
plest scalar or fermion model to intolerable degree, one 
would have to doubt that the simple constituent picture 
could suddenly become valid in QCD. If rotational sym- 
metry can be approximately obtained, one could hope 
that proceeding along the line outlined in Refs. || and 
ju^ , corrections can be refined and eventually a precise 
theory obtained. 

It was mentioned earlier in Section I that when seek- 
ing Hamiltonian dynamics for effective particles in the 
instant (or usual time) form, one would have to make 
tests of boost invariance. One could then say that the 
deviations from rotational symmetry one may observe in 
a LF approach provide a measure of how large devia- 
tions from boost covariance can be expected in models 
based on the standard form of field-theoretic constituent 
dynamics, even if they respect rotational symmetry in 
some frame of reference exactly. 



III. THE MODEL 

The gg-meson and qqg-hybrid wave functions are con- 
structed by analogy with decomposing a product of rep- 
resentations of the Lorentz group for two and three non- 
interacting particles into irreducible representations de- 
fined by their invariant mass and J quantum numbers. 
The free-particle on-shell invariant mass constraint cor- 
responds to a state with a wave function proportional to 
a (5-function in the invariant mass. In the model, mesons 
are constructed by superposing such states in the form 
of an integral with a Gaussian wave-function factor that 
depends on the invariant mass |l5f| , or, equivalently, on 
the length of the corresponding center-off-mass momenta. 
Before a theoretically sound, realistic constituent Hamil- 
tonian in QCD is identified, such approach to building 



models of the non-perturbative hadron structure is quite 
common, e.g., in the constituent models, or, with some 
variation, in QCD sum rules. 

A. qq Mesons 

The two-constituent (or qq) meson state is written as 
\p)=J2 f t 12 l P^i 1 + 2 - P)^J PC (!> 2 ) fet i4|0) , (6) 

12 ^ 

where 1 and 2 collectively represent all quantum numbers 
of the quark and anti-quark, i.e. the LF 3-momenta, p + , 
p , spins, and charges, respectively. The (free) Lorcntz- 
invariant integration measure, denoted in Eq. (^|) and 
similarly later by [12], is given by 

dk+d 2 k^ dk+d 2 k^ dx 12 d 2 k^ 2 dP+dP^ 

167r 3 fc+ 167r 3 fc+ ~ 16ir 3 a;i2(l-a;i2) 167r 3 P+ ' ^' 

where for all three components of the 3-momenta, P\ 2 = 
ki + k 2 , ki2 = (1 - xi 2 )ki - x\ 2 k 2l and x\ 2 = k\ jP^ 2 . 
Also, 5(1 + 2 - p) = 167r 3 <S(P+ - p+)S 2 {P^ - p x ). Due 
to the kinematical nature of Lorentz boost invariance 
of the LF scheme, even in the fully interacting theory, 
the coordinates P^' , x\ 2 , and k^ 2 , enable one to sepa- 
rate the relative constituent motion from the bound-state 
center-of-mass motion. In general, in the n-particle case, 
the relative-motion variables that are independent of the 
bound-state momentum p are given by (i — 1, • • ■ , n) 

x t = k+/p+, kl = ki-x iP x . (8) 

For given quantum numbers J PC , the wave function 
"J jpc (1, 2) is a product of color, flavor, spin, and orbital 
parts, and can be written as, 

(9) 

with C = for a color singlet, / = l/\/2 for an 

iso-scalar, and / = r/v2 for an iso-vector meson. Since 
iso-spin is irrelevant for our studies of Lorentz symmetry, 
further explicit consideration of the model is limited to 
iso-scalar meson states. 

The model spin-wave function x^S^l^x originates 
from the product of free Dirac spinors. In this analy- 
sis, however, we do not perform a detailed phenomeno- 
logical study for all exotic meson decay channels and 
we are not concerned about particulars of the spin-orbit 
couplings for the corresponding qq states. For example, 
J PC = o-+, gg-meson wave function could involve the 
product ui"f§v 2 as the spin factor. However, for the crit- 
ical check if the rotational symmetry can be obtained at 
all, we do not need to insist on constructing states with 
arbitrary J PC quantum numbers for the mesons in ques- 
tion (here, two, or three, non-interacting constituents). 
Instead, we choose a simplest scalar form u\v 2 for all de- 
cay products. This ansatz already introduces some spin- 
orbit structure and we can use it to check the sensitivity 
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of our numerical estimates to the presence of relativistic 
spin-orbit couplings. Note that in the calculation of de- 
cay amplitude the spins of the constituents are summed 
over and it does not matter what spin basis on chooses. 
Nevertheless, the natural choice in the LF scheme is to 
use Melosh spinors, 



u{ka) — B(k,m q )u{Qa) , 
v(k<j) = B(k, m q )v (Oct) , 
v(ka) = Cu*(ka) = i'j 2 u*(ka) , 



(10) 



and 



N p = 



2tt 3 / 4 2m q 



(17) 



with 



When we consider scalar constituents, the spin factor 
2(M.\ 2 — 4m?) is absent and the second square-root fac- 
tor in Eq. ( |13| ) i s removed, so that the normalization 
constant remains the same. 



B. Hybrid Meson 



B(k,m) 



s/k+i 



-.[A+k + +A-(m + k ± a J -)}, (11) 



which represents a LF-symmetry boost from rest to mo- 
mentum k for a free particle of mass to, using conven- 
tions A± = ±707*, u T (0a) = v^[x*\0], v T (0a) = 
y/2m q [0, £E a ], an d = —i<? 2 X<7- Adopting these con- 
ventions, one obtains 

1 



5(1,2) = 



^12(1 - X12) 



[m q (2xi2 - l)cr 3 + k^a^] ia 2 



(12) 

The orbital wave function i/> p (l, 2) = ip p (ki2) is chosen 
to be Gaussian, i.e., 



ip p (k 



12) 



IM 



12 



1 



2to, V 2(Mj 2 - Am 2 ) 



exp 



-k 2 

ft 12 

2f3 2 



(13) 

where the 3-vector kn has the _L components equal to 
k\2 and the z-component is defined by J16[, 



The construction of the hybrid meson state is analo- 
gous. 

\h) =Y^ I [123]/i + 5(l + 2+3-^)*, 7 pc(l,2,3)6 t 1 4a3l°>: 

123 

(18) 

where the color, flavor, spin, and orbital wave-function 
factors are given by 

^pc(1,2,3) = xt x Cl* Xc2 xtjhXi 2 xtS h (l,2) X s 2 
x ^(1,2,3), (19) 

with C^ 3 — t C:i /2, Ih = l/\/2, and the quark-anti-quark 
spin wave function is again of the from U1V2, so that 



S h (l,2) 



l-x 3 



Xl - X 2 3 ± j_ 

to„- (J + k„ a 

1 - x 3 q 



ia 2 . (20) 



The kinematical momentum variables satisfy Eq. (j^) 
with p + ,p ± replaced with the hybrid LF-3- momentum, 
h. Similarly, the 3-particle integration measure is given 
by, 



M 2 2 



m 2 q + kj 2 
2:12(1 - X12) 



4(fcf2 + m 2 ) , (14) 



so that k\ 2 = (x 12 — 1/2)M\2- Note that the subscript p 
in ip p {k\2) is not supposed to mean that tpp(ki 2 ) depends 
on the bound state momentum p, but rather that the 
wave function corresponds to the meson with quantum 
numbers collectively denoted by p, and the notational 
simplification should not be a source of confusion when 
we consider three mesons simultaneously. 

In terms of the 3-vector k\2, the integration measure 
over the relative motion of the two constituents, is given 
by 



16ir 3 x 12 (1 



£12) 



4d 3 fci 2 

167T 3 Xi2 



(15) 



The normalization condition (p\p') = p + 5{p — p'), im- 
plies, 



d 3 k M 



9 JVL qg 



dk+d 2 kj _ 4d 3 k q 

W 167T 3 fc+ 16TT 3 M q l 67r 3^/ m 2 + k 2 J M 2 + k 2 ' 

(21) 

where the relative momentum 3- vectors k q — (fc^~,fc 3 ) 
and kg = (kg 1 , k 3 ), are defined by 

k q = (xzk±l - Xlfe±2)/(1 - X3) , 



(22) 



and 



K = 



kf = (1 - x 3 )k£- - x 3 (ki + ki) , 
x 2 (M 2 q + k^ 2 )-(l-x 3 ) 2 (m 2 g + k^) 
2x 3 (l-x 3 )M qg 



(23) 



M q and M. qg are the quark- anti-quark and quark- anti- 
quark-gluon invariant masses, 



1 f 4d 3 k 1 



16tt 3 J M12 



2(Mi a -4m^p(l, 2)^ = 1, (16) 



M 2 = (fci + k 2 ) 2 = (1 - x 3 ) 2 m « + k « 2 , (24) 

X1X2 



6 



i 2 q + kj 

X\ 

2 + 



m 2 + k 2 



± 2 



~ 9 ' 3 -h ±2 



M 2 n + kjr 2 ml + k^ 2 



X2 
1 + 

"a 1 ,IJ g 



•i'3 



(25) 

1 - x 3 x 3 
In terms of the relative momenta, the 3-particle orbital 
wave function iph(\, 2, 3) = iph(k q , k g ) is assumed to be 



2m n 



M 



\ 



k 2 



x exp 



kq 



2(3 2 



hq 



\ 



exp 



M 2 



2m„ 



-k 2 



20? 



hg 



(26) 



Its structure follows from the requirement that it is a 
product of two Gaussian functions of the momenta k q 
and kg, with additional factors chosen so that the wave 
function normalization condition has a nonrelativistic ap- 
pearance of a six-dimensional integral over the momenta. 
In the normalization condition (h\h') = h + S(h — h'), 
fcrmion spinors introduce a factor, which is compensated 
by the factor (8fc 9 2 )~ 1//2 in the wave function. This fac- 
tor is dropped when we consider scalar constituents. So, 
in both cases of the fermion and scalar constituents the 
normalization constant is the same and equals 



since the effective gluon mass m g o is zero in the first-order 
Hqcd x ■ 

In our elementary study of the Poincare symmetry 
constraints on the constituent picture, the spin factor 
XsiSqcdO-) 2, 3)xs 2 is replaced by u\V2, which one would 
obtain in the same scheme in Yukawa theory [jlo| , and, 
using = %2ki — xik^, Sqcd(^, 2, 3) is set equal to 



5(1,2) 



1 



[m q (xi - X 2 )cr 3 + Ki 2 )cr- 



icr 



. (31) 

i.e., we remove the gluon spin effects by assuming it is a 
spin-0 particle instead. 

For convenience of normalization of the coupling con- 
stant in the test calculation, g is written as g = g\v\ 
and v\ — exp [l6m^/A 4 ] , so that for a pair produced 
at the threshold given by the constituent quark masses, 
9 fx = 9x, which is kept the same in all cases, correspond- 
ing to a x = 5a/( 47F ) = L 

In the decay, both meson states, \p) and |6) in the 
product \pb) — \p)\b), are constructed in the same way 
and may differ only trough the meson mass parameter, 
Trip or nib, and the width of the orbital wave- function tp, 
f3 p or /3b ■ Thus, the decay amplitude is a sum of two 
terms, 

A(p, b, h) = A(p, b, h) + B(p, b, h) , (32) 
that are schematically shown in Fig. [3. The result of 



N h = 



2tt 3 / 4 



I2m q 2tt 3 / 4 



8m q m g 



1 



Phq V Ph q Phg V 2m 9 + m g (3h g 



IT- ■ (27) 





C. Decay 

The constituent-gluon decay into a constituent quark- 
anti-quark pair is assumed to be driven by the Hamilto- 
nian interaction term obtained from first-order solution 
to the renormalization group equations for effective par- 
ticle Hamiltonians in quantum field theory. In QCD, the 
decay amplitude of a hybrid meson into two mesons, cal- 
culated in first order in this interaction term, has the 
form 

(pb\H IQ cDx\h) = 5(p + b-h) A(p,b,h) , (28) 
where the interaction term is given by [^| 

[123] 5(1 + 2-3) g f x (l, 2, 3) 



H 



IQCDX 



E 

123 



xlt Ca Xc 2 xlxh xISqcd(1,2,3)x 82 b\d\a 3 , (29) 

The spin-factor x1 1 Sqcd(1,2,3)Xs 2 equals u\ and 
£3 denotes the constituent-gluon polarization four-vector. 
The form factor fx is given by 

A(l, 2, 3) = fx{M\ - m 2 g0 ) = exp [-M q /X 4 ] , (30) 





FIG. 1: Hybrid meson decay amplitude into two non-exotic 
mesons p (such as 7r-meson) and b (such as 61-meson). 

summing over all quantum numbers except momenta can 
be written in the term A(j>, b, h) as 



TV 



[Sj(l, 4)S h (l, 2)5^(5, 2)S QC d(5, 4)] , (33) 



C3 

Tr 
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and in B(p, b, h) as 

Y,Tr [Cfr*alq*\ Tr[lllll h 

C3 

Tr[sl(l,4)S QC D(l, 2)4(5, 2)S h (5, 4)] . (34) 

The full expression for the amplitude including the mo- 
mentum integrals, is given in the Appendix. 

IV. NUMERICAL RESULTS 

This Section describes how the decay amplitude A de- 
pends on the angle 9 between the direction of flight of one 
of the mesons (the meson p) and the arbitrarily chosen 
z-axis, in the hybrid rest frame of reference, and how this 
dependence changes when the parameters of the model 
are varied. A true solution of QCD would not have this 
much freedom of varying parameters. However, if the 
model were found unable to satisfy the requirement of 
rotational symmetry with reasonable accuracy for any 
choice of the parameters introduced here, our effective LF 
constituent picture with only lowest Fock sectors could 
be dismissed as a good approximation to start with. 

Since the decaying hybrid and the outgoing mesons are 
constructed to be scalars, one should obtain an amplitude 
that does not depend on 9 (the LF approach respects ro- 
tational symmetry around the z-axis and all amplitudes 
we consider are automatically constant as functions of 
the azimuthal angle 4>). No dependence on 9 is always 
obtained when the sum of masses of the outgoing mesons 
is smaller by only a few MeV than the mass of the hybrid. 
In this case, both outgoing mesons move with speeds that 
are small in comparison to the speed of light. When 
at least one of the mesons becomes light and the avail- 
able mass defect forces it to move fast, the amplitude 
begins to deviate from a constant. This is illustrated in 
Fig. ^. In Fig. ||a, we illustrate what happens in the 
model with scalar constituents representing quarks, and 
Fig. ||b shows what happens in the version with quarks 
represented by fermions. The same convention of two 
versions, a) for scalar and b) for fermion constituents, is 
kept in all graphs that follow. Fig. || shows that in the 
case of one of the decay products being light, such as the 
7r-meson, the decay amplitude develops practically 100 % 
deviation from rotational symmetry. This is clearly not 
acceptable and invalidates the initial picture based on a 
conventional choice of the model parameters, given in the 
first column of Table |, labeled 2a, b. The reason is that 
when the light meson moves against the z-axis direction, 
its longitudinal momentum, 



P 



P 



(35) 
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FIG. 2: Angle dependence of the model hybrid decay am- 
plitude A for various choices of the light meson mass, given 
in MeV as indicated. The smaller is the meson mass, the 
more important become relativistic effects and the larger is 
the violation of rotational symmetry. The left plot a) is for 
scalar constituents, and the right plot b) for fermion quarks. 
All points correspond to the same standard set of parameters 
from the first column of Table 0, labeled 2a,b, and m p is varied 
from 664 MeV down to 138 MeV. See the text for details. 



becomes small. As a result of the +-momentum conser- 
vation, the quarks that make up this light meson must 
carry small fractions of the parent, hybrid momentum 



h + . This configuration is suppressed if one employs pa- 
rameters that would be naively considered as valid in a 
nonrclativistic model. In that case, the gluon decay ver- 
tex is taken directly from quantum field theory without 
the renormalization group form factor f\. 

However, varying parameters of the model one can ob- 
serve that the drastic violation of rotational symmetry is 
not always present. Since the number of mass parameters 
is 10, and varying all of them produces a huge number of 
combinations to investigate, one is faced with a problem 
of how to do a systematic test. The condition of obtain- 
ing an amplitude that does not depend on the angle 9 (an 
infinite number of conditions) is so strong that it may eas- 
ily be impossible to achieve a constant amplitude with a 
small set of parameters. It turns out, however, that one 
can obtain small variation with angle for selected sets of 
the parameters that are discussed below. 

The results are based on a systematic computer- added 
manual review of about 6000 cases and subsequent nu- 
merical minimization of the amplitude variation with re- 
spect to 9. Since the six dimensional integrals were car- 
ried out using Monte Carlo integration (using the pro- 
cedure vegas |jl7|), the accuracy of the result of integra- 
tion (standard deviation output from vegas) varied from 
a typical 0.1 - 0.3 % to 2 % in the worst cases, and 
the minimization of angle dependence (Powell's proce- 
dure pl| ) had to be carried out with a somewhat worse 
accuracy. Therefore, variations of the parameters that 
caused smaller variations of the amplitude A with 9 than 
1% could not be kept under precise control. Thus, in the 
obtained sets of parameters, displayed in Table | with 
three significant digits, the first is stable, while the sec- 
ond is highly probable and the third can be considered a 
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result of statistical fluctuation. Figs. || and |^ describe the 
two dominant effects that we found scanning the space of 
parameters, and Fig. || shows results of the final searches 
for the local minima in the angle dependence. The er- 
ror bars indicate the size of three standard deviations 
obtained in the Monte Carlo routine. 

Since the decay of a hybrid into one meson as heavy 
as b\ and another as light as tt is of key interest, the 
mass parameters = 1.9 GeV, m\, — 1.235 GeV and 
m p = 0.1375 GeV are kept constant. The other seven 
parameters were first chosen to be (all in units of GeV, 
which will be omitted in the discussion) m q = 0.3, m g — 
0.8, f3 p = /3 b = /3 hg = (3 hq = 0.4 and A = 4.0. Then a 
code calculated 3 7 cases (each for 35 different values of 9) 
with each of the 7 parameters changing through 3 values 
obtained by multiplying the initial values by 0.5, 1.0, 
and 2.0. The resulting cases were reviewed manually and 
additional calculations were performed to verify trends 
visible in the samples. This way the two generic features 
illustrated on examples in Figs. || and ^ were identified. 
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FIG. 



3: Angle dependence of the decay amplitude A with 
m p — 0.1375 MeV (corresponding to the lowest curves in 
Fig. ^) for various choices of the quark-anti-quark subsystem 
size in the hybrid state. The size is indicated in terms of 
the momentum space wave function width (3h q in units of 
GeV, as indicated (other parameters are constant and equal 
to the displayed in the second column of Table ^, labeled 3a,b). 
The spatially smaller is the pair sub-system the larger is the 
relative magnitude of the relativistic decay amplitude at large 
angles, where Fig. ^ displayed a deficiency. See the text for 
details. The left graph a) is for scalar constituents and the 
right graph b) for fermion quarks. 

Fig. U shows two cases, a) for scalar, and b) for 
fermion quarks, of what happens with the decay ampli- 
tude when one makes the quark-anti-quark pair in the 
hybrid smaller and smaller in size in position space (Phq 
grows from 1 to 5 GeV). The amplitudes are lowered 
at small angles while they remain more or less stable 
around 9 = 180°. This effect results from the fact that 
the larger domain of momenta of quarks in the hybrid 
is sampled by the same momentum-space wave functions 
of mesons, and they see less of the total probability of 
finding the quark-anti-quark pair in the hybrid, while the 



large relative momentum-tail is sampled in a similar way. 
Thus, the deficiency observed in Fig. || at large angles 
can disappear if the quark-anti-quark pair in the hybrid 
is made small in its spatial size. This is understandable 
since quarks are allowed to move with momenta so large 
that they can easily have small longitudinal fractions x 
required in the light meson, p, when it flies out against 
z-axis. However, in both cases a) and b) the amplitude 
overshoots for large 9 and a flat case is not achieved. 
This is easier to see in Fig. |I| which concerns the case 
of the lowest curves of (3} lq = 5.0 GeV from Fig. ||. Note 
that the remaining variation with angle is about thrice 
smaller in the case of fermions than in the scalar model, 
where it is still at the 100 % level. 
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FIG. 4: The result of lowering the renormalization group 
form factor width parameter, A, given in GeV as indicated, 
in the cases with m p — 0.1375 MeV and j3h q = 5 GeV (lowest 
curves in Fig. ^|). The smaller is A the more reduced is the 
amplitude at large angles. The effect confirms that lowering A 
may remove the large-angle increase in the amplitude due to 
large (5h q independently of the constituent quark spin. These 
results suggest that the small size of the quark-anti-quark pair 
in the hybrid (large f3h q ) may be correlated with the form 
factor width A in a relativistic LF-theory through rotational 
symmetry. The left graph a) is for scalar constituents and the 
right graph b) for fermion quarks. 

Fig. |J shows another tendency in the parameter space 
that is related to a reduction of the overshooting due to 
large f3hq visible in Fig. [|. Namely, the pair of quarks 
that originate from the hybrid decay via gluon dissocia- 
tion, have their relative momentum limited by the renor- 
malization group form factor f\. In a bare quantum field 
theory, A = oo and no such suppression is present. How- 
ever, when one lowers A, the gluon cannot be turned into 
a pair of very fast moving quarks through the effective 
interaction, and the amplitude is suppressed again. The 
interplay of (3hq and A turns out to be the key ingredi- 
ent of the model that at the end allows it to produce 
small deviations from a constant in the amplitude A as 
a function of 9. 

Detailed minimization of the ^-dependence leads to the 
four cases shown in Fig. [s]. The deviation from constant 
was measured using the sum of squares of deviations from 
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the average value, or maximum of the modulus of de- 
viations from the average. The values of masses, wave 
function widths, and form factor parameters that pro- 
duce the minimal angle variations, are given in Table [| 
in the last four columns. The remarkable feature is that 
the required qq-w&ve function in the hybrid has roughly 
the same width as the form factor in the effective vertex 
and this width is systematically much larger than other 
width parameters (see Table ffl). 
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ment of rotational symmetry in our relativistic LF- 
constituent picture, forces the hybrid's constituents to 
form a system made of a small-size quark-anti-quark-pair 
and a gluon, and the pair is about 5 to 10 times smaller 
than the distance between the gluon and the pair. In 
principle, this picture resembles a gluonium made of two 
gluons fH^, but with one of the gluons replaced by a 
qq-pair. The observed correlation between (3h q and A, 
suggests that the gg-pair can be seen to originate dy- 
namically from a decay of one gluon in a gluonium, a 
picture perhaps derivable from a well-defined Hqcdx- 

The amplitudes still depend on the other parameters 
than fihq and A, and those can be varied in a considerable 
ranges causing residual changes in the angle dependence. 
The most pronounced effect is that flattest amplitudes 
are systematically arrived at when quark masses take 
only one of two values, about 160-170 MeV or about 
370 MeV. Since in quantum field theory such as QCD 
one can associate the large renormalization group scale A 
with the smaller rather than the larger quark masses, it 
appears to be a remarkably consistent feature of the ef- 
fective constituent approach that the lighter quark mass 
cases lead to amplitudes that vary less in 9 than for the 
heavier quarks. 



FIG. 5: Four cases of decay amplitudes A obtained by fitting 
wave-function parameters and A so that an amplitude has 
the least possible variation as a function of the angle (see 
columns 5hs, 51s, 5hf, 51f in Table m . The two upper curves 
are obtained in the scalar model, and the two lower curves are 
obtained in the fermion model. In both models, the curves 
differ by the optimal choice of the quark mass, which is on the 
order of 160-170 MeV in the cases labeled as light, and 365-375 
MeV in the cases labeled as heavy. In all cases, f3h q ~ A ~ 4— 5 
GeV, which suggests that the size of quark-anti-quark pair is 
related to the form-factor width as if the pair in the hybrid 
state were produced through splitting of one of two gluons in 
a gluonium state of the same effective theory. See details in 
the text. 



TABLE I: Parameters in Figs. Q-|5| (in GeV). # a is associated 
with scalar particles, and b with fermions. ' hs=heavy scalars, 
ls=light scalars, hf=heavy fermions, lf=light fermions. * The 
parameter is the same as in the left neighboring column. ** 
The parameter is varying from curve to curve in the figure. 
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These results can be interpreted as that the require- 



V. CONCLUSION 

The numerical study described in this work shows that 
constraints imposed by requirements of rotational sym- 
metry on parameters of a crude LF constituent model 
of a relativistic decay of a hybrid can be satisfied even 
if the model assumes that only lowest effective-particle 
Fock-space components are included in the wave func- 
tions of the participating hadrons. Since the boost sym- 
metry is kinematical in the LF quantum field theory, and 
the model is built to be compatible with the rules of 
renormalized LF Hamiltonian dynamics for effective par- 
ticles, it is interesting that a self-consistent picture of 
the hybrid state in such a simple class of models emerges 
without major obstacles as soon as one introduces the 
right width A of the effective vertex form factor. The key 
renormalization group aspect of the constituent model is 
that the special-relativity symmetry constraints are sat- 
isfied only approximately. The accuracy is expected to 
improve when the renormalization group equations for 
effective constituent theory are solved including higher 
order terms than the first considered here, and when the 
resulting constituent dynamics is solved precisely, which 
together would be analogous to restoring rotational sym- 
metry in lattice gauge theory. 

The physical picture of a hybrid meson that one can 
now attempt to uncover by performing the detailed cal- 
culations in QCD, has the following spatial structure. 
A quark-anti-quark-pair of small size and a gluon move 
around a common center, and this system of two objects 
is not far from a gluonium. In the latter case, one has a 
second gluon, instead of the small quark-anti-quark pair. 
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Solution of the QCD Hamiltonian eigenvalue problem has 
a chance to be approximated this way only in terms of 
renormalized effective particles, for which the similarity 
renormalization group width parameter is comparable to 
the size of the wave-function parameters. The key dif- 
ficulty of the calculation will be to properly include the 
constituent-gluon spin. 
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APPENDIX A 

In the model without gluon spin, one substitutes 

SQCD(i,j) = S p (i,j) = S b (i,j) = S h (i,j) = S(i,j) , 

(Al) 

and obtains 



A(pAh) = (-l) 2 1 " f dXlid2< I 1 



3 v 7 ^ (167T 3 ) 2 J xi±{l-xu)J x 52 (l - x 52 ) x 3 
xN p N b N h r p (l, 4)^(5, 2)T(1, 2, 4, 5) [A(l, 2, 4, 5) + 5(1, 2, 4, 5)] 
16 g\ f d 3 k h f cPkn 



3 (16 



IX jdh f ^l w ^* ( fc 14) ^ ( fc 52)T(1) 2 5 4,5)[^(l ) 2,4,5) + B(l,2,4,5)] , (A2) 
7T ) J Mb J M p x 3 



where 



A(l,2,4,5) = Vft(l,2,3)/ A (A4t 5 ) , 
5(1,2,4,5) = ^ fc (5, 4, 3)/ A (M? 2 ), 



and 



T(l, 2, 4, 5) = Tr 4)5(1, 2)^(5, 2)5(5, 4)] 

with S(i,j) = UiVj, or 



(A3) 
(A4) 

(A5) 



S(i,j) 



which give T(l, 2, 4, 5) equal to 



X\X 2 X 3 Xi 



(A6) 
(A7) 



where 



m = 


[m q (xi 


X4. ) j X^k-^ 


- x\k\ 


2/2 = 


[m q (x 1 


- X 2 ),X 2 ki 


— X\k 2 


V3 = 


[m q (x 5 


- x 2 ),x 2 k$ 


- x 5 k 2 


m = 


[m q (x b 


— X4), x^k§ 


x§k^ 



(A8) 



An alternative expression is 

T(l,2,4,5) = (M^-AmDiM^-Aml) 



(M 2 15 -4m 2 g )(M 2 2 -4m 2 q ) 
(M\ 2 - 4to 2 )(M 2 5 - 4m 2 ) 



• (A9) 



In both parts of the amplitude, A and B, one has 

x p = p + /h + , Xb = b + /h + = 1 — x p . In meson p, one has 

fci4 = k p , so that Mu = M p = 2^J to 2 + fc p 2 , and the fol- 



lowing relations hold: £14 = (y to 2 , + k,^ + kp)/Mp, k^ = 

xi 4 p+, xi = xi 4 x p , kf = (1 - xi 4 )p+, x A = (1 - xu)x p , 
ki = xup ± + k p L , k^ = (l — xu)p ± — k p L . In meson 6, one 

has £52 = fcb, so that Mb = I^Jm 2 , 

analogous relations are: X52 = 

^52& + , X 5 = X 52 X b , fc+ 



+ k b = A^52, and the 



(Vm 2 + k 2 + k 3 b )/M b , 
(1 - x 52 )b+, x 2 



(1 - x 52 )x bl k^ = x^b 1 - + k^, k^ = (1 - x 52 )6 ± - k^- 



b ■ 



Evaluating the quarks' invariant masses in the hybrid and 
decay vertex, one then obtains 



M 2 2 = 



{xi +x 2 ) 
(kt + ki) 2 



kt 



xi 



+ 



K 2 



Xi 



(A10) 



and 



M\± = (x 5 + x^ 



H 2 



ki 2 



in,, 



x 5 



In the part A of the decay amplitude, one has 



k 2 = M\ 2 /A - to 2 



(All) 



(A12) 
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and 



and 



k: = 



where 



MI23 ~ (M12 + m g ) M\ 23 - {M12 ~ m g ) 



2:i 



k^ 2 + m 2 k,2 2 + m 2 



(A13) 



.Tl 



X 2 



(to + to) X2 +m 2 a , . „ 
= ^ - (J* + fe 2 + to) ±2 



1 - £l - £2 



Similarly, in the part £?, one has 



(A14) 



k 2 = 



where 



M 2 5i3 - (M 5 A + m g 



{M 5 4 - m g Y 



4ML 



(A16) 



A4 



5 VI 



k^ 2 + m 2 k^ 2 + m 2 
(k 5 + k i ) X2 +m 2 g 

-(to + to + to) 



1 — X5 — X4 



(A17) 



k 2 =M\j±-m\ , 



(A15) 
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